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Optimization of tilt-rotor systems requires the consideration of performance at multiple 
design points. In the current study, an adjoint-based optimization of a tilt-rotor blade is con- 
sidered. The optimization seeks to simultaneously maximize the rotorcraft figure of merit 
in hover and the propulsive efficiency in airplane-mode for a tilt-rotor system. The design 
is subject to minimum thrust constraints imposed at each design point. The rotor flow- 
fields at each design point are cast as steady-state problems in a noninertial reference frame. 
Geometric design variables used in the study to control blade shape include: thickness, cam- 
ber, twist, and taper represented by as many as 123 separate design variables. Performance 
weighting of each operational mode is considered in the formulation of the composite objec- 
tive function, and a build up of increasing geometric degrees of freedom is used to isolate the 
impact of selected design variables. In all cases considered, the resulting designs successfully 
increase both the hover figure of merit and the airplane-mode propulsive efficiency for a rotor 
designed with classical techniques. 


Nomenclature 


c 

Aerodynamic coefficient 

S 

Control volume surface area 

c 

Chord 

T 

Temperature 

D 

Vector of design variables 

t 

Time 

E 

Total energy per unit volume 

U 

Velocity 

f 

Objective function 

U, V, w 

Cartesian components of velocity 

F, 

Inviscid flux tensor 

Ui 

Cartesian directional displacements 

FM 

Rotorcraft figure of merit 

V 

Volume of control volume 

F 

1 V 

Viscous flux tensor 

X 

Vector of mesh coordinates 

1 

Identity tensor 

Xi 

Cartesian coordinate directions 

K 

Elasticity coefficient matrix 



k 

Thermal conductivity 

Subscripts 

N 

n 

P 

P 

Q 

R 

r 

R 

r 

S 

Number of blades 

Outward-pointing unit normal vector 
Pressure 

User-defined exponent 
Vector of conserved variables 
Spatial residual vector 
Position vector 
Blade radius 
Radial station 
Source term vector 

cruise 

e 

f 

hover 

i,j,k 

m 

mp 

Q 

surf 

T 

Cruise condition 
Effective 
Flow field 
Hover condition 
Indices 
Mesh 

Multi-point composite 
Torque 

Surface quantity 
Thrust 
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tip Rotor tip 

Symbols 

0 () Order of Magnitude 

a User-defined weighting factor 

A Difference 

e Strain tensor 

77 Propulsive efficiency 

O Blade collective pitch setting 

A Adjoint variable vector 

X, p Linear elasticity Lame constants 


v Poisson’s ratio 

p Density 

a Rotor solidity 

r Viscous stress tensor 

$7 Angular velocity vector 

ui User-defined weighting factor 
00 Freestream quantity 

Superscripts 

T Vector transpose 

* User-defined target 


I. Introduction 

A pplication of high-fidelity computational fluid dynamics (CFD) has become commonplace in the fixed- 
wing aerospace community. Software packages that solve the Euler equations and Reynolds-averaged 
Navier-Stokes equations on both structured and unstructured meshes are now used routinely by aerodynami- 
cists in the analysis and design of new configurations. Moreover, as algorithms and computer hardware have 
continued to mature, the coupling of formal design optimization techniques and sufficiently accurate physical 
models has become a viable solution for the design of many large scale aerospace configurations. 

The application of high-fidelity CFD tools to the analysis and design of full rotorcraft configurations 
is considerably challenging. Such flowfields are inherently unsteady, frequently involve fluid velocities that 
range from quiescent to transonic flow, and typically require the simulation of complex aerodynamic and 
aerostructural interactions between dynamic vehicle components. Recent literature suggests that the use of 
high-fidelity CFD methods in this regime is growing, but the computational cost that is required to capture 
the necessary spatial and temporal scales of a typical rotorcraft flowfield remains considerable. 1-9 

In the field of gradient-based design, adjoint methods are known to provide an extremely efficient means 
for computing sensitivity information. The cost of such methods is equivalent to the expense that is associated 
with solving the analysis problem and is independent of the number of design variables. Adjoint methods can 
also be used to perform mathematically rigorous mesh adaptation and error estimation. Significant success 
has been reported for the application of these techniques to steady problems. 10-14 

In general, optimization and mesh adaptation for large-scale unsteady flows based on adjoint methods 
require a time-dependent implementation of the equations. Considerable effort by a number of research groups 
is being focused in this area, and examples of the use of such approaches have recently emerged. 15 ~ 17 Despite 
the algorithmic efficiency, however, the computational cost of these general time-dependent approaches can be 
considerable. 

An adjoint-based design capability was previously demonstrated for rotor configurations for which the 
analysis problem may be cast as a steady problem in a noninertial reference frame. 18 This approach permitted 
the use of an existing steady-state adjoint formulation with minor modifications to perform sensitivity analyses. 
The resulting formulation was valid for isolated rigid rotors in hover or where the freestream velocity is aligned 
with the axis of rotation. 

The goal of the current work is to apply the non-inertial capability 18 to the multi-point design of a tilt-rotor 
blade. The design will simultaneously optimize aerodynamic performance of the rotor at both hover and cruise 
conditions. For hover, the objective is to maximize the rotor figure of merit while maintaining a minimum 
thrust coefficient. For cruise, the objective is to maximize the cruise performance efficiency while maintaining 
a minimum thrust coefficient. To better understand the contributions of the chosen geometric design variables, 
the design follows an orderly build up of increasing geometric degrees of freedom. 

II. Flow Equations 

The governing equations for the flowfield are the compressible, perfect gas Reynolds-averaged Navier-Stokes 
equations written in a reference frame that is rotating with a constant angular velocity, L7: 


d(QV) 

dt 



F„) • n dS = S 


(1) 
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where Q is the vector of volume-averaged conserved variables, Q = \p, pu , pv , pw , E] T , n is an outward-pointing 
unit normal, and V is the control volume bounded by the surface dV. The inviscid and viscous flux tensors 
are given by 


F- 


p (u — i"2 x r) 
pu T (u - U x r) + PI 
(u — ft x r) (E + P) + (fi x r) P 


( 2 ) 


and 
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0 

T 

u ■ r — fcVT 


The source term S represents a Coriolis effect that results from the rotating frame of reference: 


( 3 ) 


S = 


0 

— p ($1 x u) 
0 


( 4 ) 


Here, u is the absolute velocity vector, u = [u, v, w) T , r is the position vector relative to the axis of rotation, and 
r is the viscous stress tensor. The equations are closed with the perfect gas equation of state and an appropriate 
turbulence model for the eddy viscosity. For rotorcraft simulations, the formulation that is described here is 
applicable to rigid tilt-rotor geometries in either a hover, ascending/descending, or airplane mode cruise flight 
condition, where the freestream velocity vector is parallel to the angular velocity vector, ft. 

The FUN3D flow solver that is used in the current work 19-22 can be used to perform aerodynamic simula- 
tions across the speed range, and an extensive list of options and solution mechanisms is available for spatial 
and temporal discretizations on general static or dynamic mixed-element unstructured meshes that may or 
may not contain overset mesh topologies. 

In the current study, the spatial discretization uses a finite-volume approach in which the dependent 
variables are stored at the vertices of single-block tetrahedral meshes. Inviscid fluxes at cell interfaces are 
computed using the upwind scheme of Roe, 23 and viscous fluxes are formed using an approach that is equivalent 
to a central difference Galerkin procedure. The eddy viscosity is modeled using the one-equation approach 
of Spalart and Allmaras 24 with the source term modification proposed by Dacles-Mariani et al. 25 For the 
steady-state flows (relative to the noninertial reference frame) that are described in this study, temporal 
discretization is performed using a backward-Euler scheme with local time stepping. Scalable parallelization 
is achieved through domain decomposition and message-passing communication. 

An approximate solution of the linear system of equations that is formed within each time step is obtained 
through several iterations of a multicolor Gauss-Seidel point-iterative scheme. The turbulence model is inte- 
grated all the way to the wall without the use of wall functions. The turbulence model is solved separately 
from the mean flow equations at each time step with a time integration and a linear system solution scheme 
that is identical to that employed for the mean flow equations. 


III. Mesh Equations 

To deform the interior of the computational mesh as the surface mesh evolves during a shape-optimization 
procedure, the mesh is assumed to obey the linear elasticity equations of solid mechanics. These relations can 
be written as 



I-ndS- 


2 p£ ■ n dS = 0 


dv 


( 5 ) 
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where 


1 ( dui duj \ 

2 \ dxj dxi J 


( 6 ) 


is the strain tensor, iq is the displacement vector in each of the Cartesian coordinate directions, Xi, and A and 
/r are material properties of the elastic medium. The quantities A and fj, are related to Young’s modulus, E, 
and Poisson’s ratio, v, through the following: 


A = 


uE 

(1 + v) (1 — 2u) 


( 7 ) 


and 


E 

2(l + i/) 


( 8 ) 


The system is closed with the specification of two of the four parameters A, fi, E 1 and v. In the current 
implementation, E is taken as inversely proportional to the distance from the nearest solid boundary, while 
Poisson’s ratio is taken uniformly as zero. This approach forces all cells that are near boundaries to move in a 
nearly rigid fashion, while cells that are far from the boundaries are allowed to deform more freely. The system 
of equations is solved using GMRES 26 with either a point-implicit or ILU(O) preconditioning technique. 22,27 


IV. Discrete Adjoint Equations 

To derive the discrete adjoint equations, a compact notation is introduced for the governing equations that 
are outlined above. The spatial residual vector R of Eq. (1) is defined as 


R = <£ (Fi - F„) • n dS - S (9) 

JdV 

Furthermore, the linear system of equations given by Eq. (5) can be written as 


KX — X. sur j 


(10) 


where K is the elasticity coefficient matrix that results from the discretization of Eq. (5) , X is the vector of 
the mesh-point coordinates, and X sur f is the vector of known surface mesh-point coordinates, complemented 
by zeros for all interior coordinates. With the approach that was taken by Nielsen, 11 a Lagrangian function 
can be defined as follows: 


L (D, Q,X, A/, A m ) 


/ (D, Q, X) + AjR (D, Q, X) + \ T m (KX - X surf ) 


( 11 ) 


where D represents a vector of design variables, / is an objective function, and A f and A m are the adjoint 
variables that multiply the residuals of the flow and the mesh equations. In this manner, the governing 
equations may be viewed as constraints. 

Differentiating Eq. (11) with respect to D and equating the <9Q/<9D and 5X/9D coefficients to zero yields 
the discrete adjoint equations for the flowfield and mesh, respectively: 

<9R 

5Q 
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and 


K 1 A„, = - 



dR 

SX 


(13) 


The remainder of the terms in the linearized Lagrangian can be grouped to form an expression for the final 
sensitivity vector: 


dL _ df . r _ * t 
dD <9D + f <9D 

Equations (12) and (13) provide an efficient means for determining discretely consistent sensitivity information. 
The expense that is associated with solving these equations is independent of D and is similar to that of the 
governing equations. After the solutions for Af and A m have been determined, then the desired sensitivities 
may be calculated using Eq. (14), for which the computational cost is negligible. 

A discrete adjoint implementation has been developed 11, 17,21,27,28 for the flow solution method that is 
described above. The flowfield adjoint equations are solved in an exact dual manner, which guarantees an 
asymptotic convergence rate that is identical to the primal problem and costate variables that are discretely 
adjoint at every iteration of the solution process. The grid adjoint equations are solved using GMRES in a 
manner that is identical to the method used for Eq. (5) . To accommodate the noninertial reference frame used 
in the current study, minor modifications have been made to include the effects of the mesh speeds and the 
Coriolis terms. 


d^-surf 

<9D 


(14) 


V. Design Methodology 


A. Design Variables 

The implementation that is described by Nielsen 11 is sufficiently general such that the user is able to employ 
a geometric parameterization scheme of choice, provided that the associated linearizations required by the 
adjoint method described above are also available. For the current study, a mesh parameterization scheme 
tailored for aerodynamic shape optimization is used. 29 This approach can be used to define very general shape 
parameterizations of existing grids using a set of aircraft-centric design variables, such as camber, thickness, 
shear, twist, and planform parameters at various locations on the geometry. The user also has the freedom 
to directly associate two or more design variables to create more general parameters. In the current work, 
this option is used to link several piecewise twist variables across the span of a rotor blade to create a single 
twist variable that is used to prescribe the blade collective pitch setting, 0. The option is used similarly to 
link planform variables to control blade taper. In the event that multiple bodies of the same shape are to be 
designed, as in the case of rotor geometries, the implementation allows a single set of design variables to be 
used to simultaneously define all bodies as a group. In this manner, all bodies in the group have the same 
shape at a given design cycle throughout the course of the design. 

B. Objective and Constraint Functions 

The implementation of Nielsen 11 permits multiple objective functions, fi, and explicit constraints, Cj , of the 
following form, each containing a summation of n, and mj individual components, respectively: 


Ili 

/, = 5> fc (c fc - c fc *r (15) 

k = 1 

and 

mj 

c,- = (C fe - C fc *r ( 16 ) 

k=l 
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Here, Uk represents a user-defined weighting factor, Ck is an aerodynamic coefficient such as total drag or 
the pressure or viscous contributions to such quantities, and pk is a user-defined exponent. The * superscript 
indicates a user-defined target value of Ck- Furthermore, the user may specify the boundaries in the grid to 
which each component function applies. 

C. Design Points and Optimization Strategies 

The current implementation supports an arbitrary number of user-specified design points at which objective 
and constraint functions may be posed. Each design point may be defined by a variation of basic flowfield 
quantities such as the Mach number, or a more general characteristic such as a specific computational grid 
that is appropriate for each individual design point. In the current study, the two design points use their own 
distinct baseline mesh, one for hover and the other for the cruise condition. While it would have been possible 
to elastically deform a single mesh to accommodate the changes in 0 for the two conditions, the individual 
meshes were generated to insure the highest quality baselines about which subsequent design driven shape 
changes would be incorporated through elastic deformation. 

To perform multi-point optimization, a composite objective function, / mp , can be defined based on the 
individual objective functions, /* , posed at each design point. In this work, the composite objective function, 
Imp, defines a linear combination of the objective functions for the hover and cruise conditions: 


fmp — Oh, over fh over H "~ 1 


l cruiseJ cruise 


(17) 


where a is a constant weighting factor for each point that is applied to the corresponding objective at the 
point. 

The multi-point approach that is used here is commonly used to combine point solutions to multi-objective 
optimization problems into a single scalar objective. The difficulty is that out of the range of many possible 
solutions only one is obtained by setting some parameters heuristically and externally, for example the weights 
of the composite scalar objective. The investigation of more sophisticated optimization strategies is out of 
scope of the current work. In principle the simple strategy adopted here is sufficient for the problem of interest 
and the investigation of more sophisticated optimization strategies is relegated to future work. 

SNOPT 30 ’ 31 is used to solve the optimization problem defined above subject to explicit constraints specified 
at each design point. The optimization algorithm is allowed to perform a maximum of 20 design cycles with a 
maximum of 30 function evaluations per cycle, and considers the design converged and exits if it believes the 
current objective function matches the value at the optimal solution to four significant digits. Constraints are 
considered satisfied if their values do not exceed the specified bounds by 0.1% of the bound value. The design 
at the initial choice of D is not required to satisfy the constraints; if needed, the optimizer attempts to locate 
a feasible starting point on its own. 

D. Design Case 

The subject of the current study is the three-bladed Tilt Rotor Aeroacoustics Model (TRAM) 32,33 shown in 
Figure 1. The TRAM represents an optimized system designed with traditional tools. It was selected for this 
study because the model readily exists and is an adequate surrogate for the more sensitive V-22. The original 
quarter scale TRAM model was resized to match full scale flight conditions. The multi-point optimization 
considers both hover and airplane mode cruise conditions similar to those found in open literature for the V-22. 
In hover, the tip Mach number is 0.707, the freestream temperature is 519° Rankine, and the Reynolds number 
is 9.2 million, based on the blade tip chord of 22 inches. The cruise condition was modeled at an altitude of 
3,000 feet, a tip Mach number of 0.650, a freestream temperature of 551.5° Rankine, and a Reynolds number 
of 6.1 million (again based on the blade tip chord). 

A separate baseline mesh was generated for each of the two conditions of the design study. The mesh for 
the hover evaluation contained 8,105,037 nodes and 47,884,674 tetrahedral elements and was designed for a 
nominal collective pitch setting 0 of 14°. The cruise condition mesh contained 8,106,096 nodes and 47,891,211 
tetrahedral elements at a nominal 0 of 41 °. The same relative spacing constraints were used to generate both 
meshes as indicated by the similar sizes of the meshes. The surface grid for one of the blades is shown in 
Figure 1(b). All of the grids have been generated with VGRID. 34 The blade trailing edges are blunt. 
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(a) Rotor surface geometry (b) Typical blade surface mesh used for design com- 

putations 

Figure 1. Tilt Rotor Aeroacoustics Model (TRAM) 


A geometric parameterization 29 has been developed for the baseline blade geometry, as shown in Figure 2. 
The approach yields a total of 123 active design variables, including 50 variables to control the blade thickness, 
60 variables to control the blade camber, 10 to control the blade twist, and 1 to control the blade taper. 
The final two variables were defined to address the collective pitch angles, 0, for each design point. These 
were formed by linking together all of the blade twist variables at a given design point into a single design 
variable. A separate design variable was required for each design point as control of 0 must be independent 
at each point. The taper design variable was also formed by linking together the two in-plane components 
of the planform variables at the blade tip. However, in this linking, the components were linked such that 
the movement of the aft planform variable was three times the magnitude of the forward variable and in the 
opposite direction. This was defined such that the |-chord line would remain invariant to changes in taper. 
To further facilitate the invariance, the |-chorcl line of the blade was aligned to the parameterization reference 
frame through a series of rigid body transformations prior to parameterization. 

Bounds on the design variables have been initially chosen with the intent to prevent nonphysical surface 
shapes; only thickness changes that increased the section thickness were allowed to eliminate thinning of the 
rotor blades. The parameterization also allows for radial blade shearing deformations, however these variables 
are not active design variables in the current study. All blades at a given design point are grouped together 
such that a perturbation of the design variables produces the same blade shape for the rotor system. Other 
than the collective pitch angle, 0, a common parameterization is used between design points. This maintains 
the same blade shape across the design. 


• Camber and Thickness 
° Camber 


T|=0.15 r|=0.25 r|=0.35 r|=0.45 r|=0.55 r|={).65 r|=0.75 r|=0.85 r|=0.95 



Figure 2. Design variable and radial blade locations. 


The design problem definition calls for explicit constraints on the thrust coefficient, Ct , at each condition. 
Prior to the multi-point design, a simple preliminary design optimization was performed for each condition 
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to determine the 0 required to trim the rotor to the target Ct- The trimming was performed as two sep- 
arate single-point unconstrained optimizations where 0 was the only design variable and the respective Ct 
represented a single objective for each point. The resulting changes to 0 were modest, as will be shown 
below. Again, it is not required that the constraints be satisfied to begin the optimization, however, this did 
demonstrate the ability of the optimizer to trim the rotor to a target thrust. 

For the hover condition, a single objective is used, where u = 1, p = 2, and C is defined as the square of 
the commonly used rotorcraft figure of merit, FM , which is composed of the rotor thrust, Ct, and torque, 
Cq, coefficients: 


C = FM 2 


C f. 
2C| 


(18) 


The square has been introduced to avoid the appearance of a square root of a negative thrust value in the 
linearized form of the objective function. The value of C* , of Eq. (15), is chosen to be 2, which is considerably 
larger than both the baseline value and the theoretical maximum value of 1. By defining C* and p in this 
way, minimization of the objective will minimize the difference between the actual and target coefficients in 
Eq. (15) effectively maximizing the value of FM. The explicit thrust constraint was used to set a lower bound 
on Ct of 0.0158 with lu chosen so that the constraint was of 0( 1). 

The cruise condition also used a single objective function with u> = 1, p = 2, and C defined as the airplane 
mode propulsive efficiency rp. 


_ _ C T Uoo 

~ V ~ CqU Up 


(19) 


Again, the value of C* is chosen to be 2, which is considerably larger than both the baseline value and the 
theoretical maximum value of p = 1. Minimization of the objective will again maximize p. The explicit thrust 
constraint placed a lower bound on Ct of 0.0027 where u> was again chosen so that the constraint was of 0(1). 

All computations were performed using the NASA Advanced Supercomputing (NAS) high performance 
computing resources. A design cycle requires at a minimum a single function and gradient evaluation for the 
current value of D. A function evaluation in this context consists of an evaluation of the surface parame- 
terization for each blade, a solution of Eq. (5) to deform the interior of the mesh according to the current 
surface grid, and a solution of the flow equations, Eq. (1). Using the adjoint approach that is outlined above, 
a gradient evaluation requires a solution of the flowfield adjoint equations, Eq. (12); a solution of the mesh 
adjoint equations, Eq. (13); an evaluation of the linearized surface parameterization for each blade; and, finally, 
an evaluation of the gradient expression given by Eq. (14). The convergence criteria that are used for each 
of the solvers has a direct impact on the efficiency. Note, the time that is required to solve Eqs. (1) and (12) 
tends to decrease toward the end of an optimization as the design converges and the solution restarts become 
more effective. 

E. Design Considerations 

In tilt-rotor design, rotor efficiency must often be a trade off between airframe and wing efficiency while 
aerodynamic efficiency is traded against weight and mission efficiency (e.g. total fuel burn). This leads to 
a paradox: maximizing aerodynamic efficiency of an isolated rotor does not guarantee maximum operating 
efficiency, and may even lead to a poorly designed aircraft. 

In hover, tilt-rotor thrust must match weight plus download. In cruise, thrust must match drag, including 
induced drag from the airframe and wing. Consider a cruise velocity determined by mission requirements and 
a wing design optimized for that velocity. If the mission is fixed and the aircraft is not resized, then neither 
hover nor cruise thrust should change. Any improvement provided by aerodynamic optimization should be 
manifest as a reduction in power (equivalently represented 35 by Cq in Eqs. (18) and (19)), not an increase in 
thrust or airspeed. 

Furthermore, while in hover most of the rotor power results from induced power. However, in cruise the 
power is primarily dependent on profile power, which varies little with thrust when a well designed rotor is 
operating at its design condition. Therefore, in high-speed cruise, Ct can be increased with negligible increase 
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in profile power, resulting in a large increase in 77. However, the increased thrust cannot be effectively utilized 
by an aircraft operating at fixed velocity given a fixed value of airframe drag. Alternatively, increasing the 
velocity at fixed power may require a heavier wing to satisfy aeroelastic stability margins. As a result, increases 
in drag or weight that result from higher cruise Ct may result in worse overall mission performance than a 
configuration with lower Ct and lower rj. Hence the above paradox that maximizing rotor efficiency may not 
maximize operating efficiency. 

The argument may be inverted when considering the trade-offs between cruise and hover efficiency. The 
aircraft may benefit from a smaller rotor with higher cruise Ct, which would give a lower FM but higher 77. 
Such a result should not be assumed in advance, however. The optimum design may require a larger rotor with 
higher FM and lower 77 . In either case, rotor performance trends for both hover and cruise must be considered 
based on a consistent reference rotor size, which is determined by the thrust needed for hover. 

The above discussion is predicated on the assumption that one sizes the aircraft to some combination of 
hover and cruise efficiency. But there are other design constraints not mentioned, notably engine-out and 
maneuver requirements. Design for engine-out conditions is implicitly included in a design to maximum FM. 
Design for engine-out requires the maximum possible thrust at a strict power limit, which also maximizes FM 
by definition. Design for maneuvers, however, generally requires maximization of available thrust, which will 
not occur at maximum FM. Worse, critical maneuvers, particularly low-speed turns, do not have axial flow at 
the rotor, even for tilt-rotors or tilt-wings. In turning flight, the total lift required is greater than the weight. 
For tilt-rotors, the worst-case turning condition will occur at low speed with the rotor in edgewise flight or 
tilted slightly forward, and with the wing carrying only part of the total lift. 

As such, aerodynamic optimization based on thrust, as distinct from maximum efficiency, is nevertheless 
required. If greater thrust can be achieved at a given rotor size, then the aircraft can be resized with a smaller 
rotor. A smaller rotor will have both lower profile power and higher Ct in cruise, which would automatically 
provide higher 77 without compromising airframe aerodynamics. A sizing code can be used to determine the 
proper trade-off between maximizing thrust in maneuvers and maximizing FM in hover with aerodynamic 
optimization via CFD utilizing a consistent rotor size. For acceptable CPU time, and to minimize problems 
with convergence, an axial-flow surrogate for a turning condition can be used whereby a constraint is added 
that requires some minimum value of thrust at nominal operating conditions, where the specified Ct is greater 
than that at which FM is maximized. 

Thrust can be increased by increasing blade chord, but this does not increase FM and decreases 77 because 
profile drag is higher. To properly separate the effects of blade area from the effects of optimized airfoils, 
twist, and planform, rotor performance coefficients should be scaled by the ratio of blade area to disk area, 
or solidity, a. Local blade section velocity, dynamic pressure and Mach number all increase with radius, so 
designing rotors to thrust-weighted solidity 36 is common: 


where 


_ Nc e 
nR 

(20) 

Jo* cr2dr 

Jo* r2dr 

(21) 


At the time of this writing, the full linearization of the blade loading (Ct/ct), namely the linearization of u, 
was not available. Therefore, it was not possible to use blade loading as an explicit constraint in this design. 
As an alternative, we imposed a minimum constraint on Ct as mentioned above for each mode of flight. 


F. Design Breakdown 

To better understand the contributions of airfoil shape, twist, and planform toward optimized rotor design, 
the design was broken into four distinct steps referenced A-D. Step A sought to optimize airfoil camber 
and thickness distributions at radial increments along the blade to maximize hover FM and cruise rj while 
keeping airfoil twist and planform fixed. A hover FM to cruise performance 77 weighting of 2:1 was selected 
by setting the weights of the composite objective function of Eq. (17) as ahover = 2.0 and a cru i se = 1.0. Step 
B added airfoil twist as an additional design variable to Step A while maintaining the same 2:1 hover-to-cruise 
performance weighting. Step C built upon previous steps to add blade taper. Again the 2:1 ratio of FM to 77 
was maintained. Finally, Step D repeated Step C with the hover-to-cruise weighting reversed from 2:1 to 1:2, 
Q-hover = 1-0 and a. cru i se = 2.0. The design steps are summarized in Table 1. 
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Table 1. Breakdown of design steps 


Step 

Camber 

Thickness 

Twist 

Taper 

^ hover 

^ cruise 

A 

/ 

/ 



2.0 

1.0 

B 

/ 

/ 

/ 


2.0 

1.0 

C 

/ 

/ 

/ 

/ 

2.0 

1.0 

D 

/ 

/ 

/ 

/ 

1.0 

2.0 


VI. Results 

The multi-point design objective was to maximize both the FM in hover and 77 in cruise. The optimization 
was simultaneously required to maintain or exceed the baseline thrust coefficients for each operational mode. 
In addition to the shape variables active for each design step, a decoupled collective pitch design variable was 
defined at each design point resulting in two additional design variables per step. 

All optimizations were executed on 2,048 Intel® Sandybridge processing cores of the NAS Pleiades system. 
Note that the wall time is not a linear function of number of function evaluations due to many complicating 
factors (each function evaluation may require a different number of iterations to converge, a design cycle may 
require multiple function/gradient evaluations, etc.). Step A activated 110 shape variables to control airfoil 
camber and thickness. With the addition of a 0 variable for each of the two design points, the optimization of 
112 variables required 17 function evaluations and roughly 20 hours of wall time. The addition of blade twist 
in step B resulted in 120 shape variables. A total of 48 function evaluations were required to optimize the 122 
variables, and roughly 49 hours of wall time. 

Both Steps C and D employed all 121 active shape variables to optimize airfoil camber, thickness, twist, and 
blade taper. Addition of the decoupled 0 design variables defined for each design point resulted in 123 design 
variables in total. In these remaining steps, the composite design objectives were subject to two inversely 
proportional performance weightings. Step C required 29 function evaluations and roughly 53 hours of wall 
time. Step D required 51 function evaluations and roughly 43 hours of wall time. 

A. Hover Design Point 

Tables 2 and 3 detail the output from each step in the design buildup for the hover point. For each step in 
the tables, the optimized collective pitch angle, 0, is given along with the values of the thrust constraint, Ct, 
blade loading Ct/<t, and figure of merit, FM. Also shown is the difference in the figure of merit, A FM, and 
its percent change from the baseline value, %Change. 


Table 2. Figure of Merit function (Constant Ct) 


Step 

® hover 

C T 

Ct 

a 

FM 

A FM 

%Change 

Baseline 

13.91° 

0.0158 

0.1505 

0.7564 

- 

- 

A 

13.52° 

0.0159 

0.1514 

0.7685 

0.0121 

1.600 

B 

14.60° 

0.0158 

0.1505 

0.7640 

0.0076 

1.005 

C 

13.83° 

0.0158 

0.1365 

0.7707 

0.0143 

1.891 

D 

14.22° 

0.0158 

0.1386 

0.7603 

0.0039 

0.516 


Table 3. Figure of Merit function ( Ct increase allowed) 


Step 

® hover 

C T 

Ct 

a 

FM 

A FM 

%Change 

Baseline 

13.91° 

0.0158 

0.1505 

0.7564 

- 

- 

A 

13.47° 

0.0169 

0.1610 

0.7714 

0.0150 

1.983 

B 

11.61° 

0.0160 

0.1524 

0.7828 

0.0264 

3.490 

C 

12.85° 

0.0175 

0.1563 

0.7840 

0.0276 

3.649 

D 

12.17° 

0.0168 

0.1556 

0.7819 

0.0255 

3.371 


Table 2 relates the cases where an equality constraint was imposed to hold Ct at or near its baseline value. 
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Table 3 shows data from the optimization where Ct was allowed to vary above the baseline value. Recall 
Section V, E explained why arbitrary increases in Ct may not result in an good design. However, in the 
absence of a blade loading constraint, the data of Table 3 show that increases in Ct, such as would accompany 
maintaining Ct/ct for steps C and D, result in further increases in FM. 

Also note that the trimmed baseline value of 0 is quite close to the nominal value of 0 = 14° used to 
develop the mesh. Each design step showed a modest increase in the FM with the greatest increase shown for 
Step C. Considering Table 3, it is expected that an equality constraint on Ct/ct holding to the baseline value 
of 0.1505 would produce better than a 2.0% increase in FM. 

B. Cruise Design Point 

Tables 4 and 5 represent the results for the design buildup at the airplane-mode cruise point. Here the 
optimized collective pitch angle, 0, is given along with the values of the thrust constraint, Ct, blade loading, 
Ct/ct, and propulsive efficiency, rj. Also shown is the difference in the propulsive efficiency, At/, and its percent 
change from the baseline value, %Change. 

Similar to the hover data, Table 4 represents the cases where an equality constraint was imposed to hold 
Ct at or near its baseline value of 0.0027. Table 5 shows data from the optimization where Cr was allowed 
to vary above the baseline value. The impact of increased Ct on 77 is clearly shown in Table 5. While not 
all aircraft can exploit the increased 77 shown in Table 5, the results illustrate the potential gains in cruise 
efficiency where the aircraft design space allows increases in Ct- 

Table 4. Propulsive efficiency (Constant Ct) 


Step 

® cruise 

Ct 

Ct 

a 

1 

Arj 

%Change 

Baseline 

41.37° 

0.0027 

0.0257 

0.6738 

- 

- 

A 

40.98° 

0.0028 

0.0267 

0.6816 

0.0078 

1.158 

B 

41.78° 

0.0027 

0.0257 

0.6971 

0.0233 

3.458 

C 

41.71° 

0.0027 

0.0233 

0.7023 

0.0285 

4.230 

D 

42.05° 

0.0027 

0.0234 

0.7008 

0.0270 

4.007 


Table 5. Propulsive efficiency (Ct increase allowed) 


Step 

® cruise 

Ct 

Ct 

a 

1 

Ag 

%Change 

Baseline 

41.37° 

0.0027 

0.0257 

0.6738 

- 

- 

A 1 

42.50° 

0.0076 

0.0724 

0.8243 

0.1505 

22.34 

B 1 

42.50° 

0.0098 

0.0933 

0.8491 

0.1753 

26.02 

C 1 

42.50° 

0.0089 

0.0795 

0.8469 

0.1731 

25.69 

D 1 

42.50° 

0.0098 

0.0907 

0.8506 

0.1768 

26.24 


1 © at Upper bound 


Here too, the trimmed baseline value of 0 is near the nominal 0 = 41 ° used to construct the reference 
mesh. Each design step showed an increase in the propulsive efficiency with the greatest increases shown for 
Steps C and D where the planform was increased. It is reasonable to assume from this data that increases in 
77 in excess of 4% for Steps C and D can be expected as Ct/ct is allowed to approach the baseline value of 
0.0257%. 

Note that for each design step when Ct was constrained to be greater than the baseline value, the collective 
pitch angle 0 was increased to its upper bound for the design of 42.5 °. This accounts for much of the increase 
in Ct and further points to the need for a true Ct/ct. 

C. Design Results 

Both the hover and cruise design points show modest increases in both the FM and rj with increasing geometric 
degrees of freedom. The weighting used in Step C increased the hover FM a modest 1.89% over the baseline 
value for fixed Ct, while the cruise performance increased 4.23% over the baseline. When inverting the 
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weighting for Step D to 1:2, the hover figure of merit improvement dropped to 0.52% while the cruise r ; 
maintained improvement at 4.0%. 

Figures 3-6 represent airfoil sections from r/R = 0.20 to r/R = 1.0 for each of the design Steps A-D 
respectively. Note that perturbations to the blade collective pitch 0 have been removed to allow comparison 
to the baseline shapes that are represented by the solid black lines. Also note that the vertical scale has been 
exaggerated to emphasize the shape differences. The airfoil shapes representing the design configuration are 
represented with dotted red lines in each figure. 

In Figure 3, modest shape changes are observed for Step A with the most notable changes occurring 
outboard of the r/R = 0.50 section. Increases in camber are observed near the trailing edge for sections 
r/R = 0.50 through 0.95. Thickness changes are most notable for sections r/R = 0.60 to 0.95. Figure 4 
shows the same airfoil sections for Step B. The shape changes resulting from the addition of the blade twist 
variable are more readily apparent and are now observed over the entire blade span from r/R = 0.20 to 1.0. 
In addition, the rotor blade appears to exhibit some washout near the tip from r/R = 0.90 to 1.0. 

The optimized airfoil sections resulting from Step C are shown in Figure 5. The addition of blade taper 
as a design variable is reflected in the increased chord of the outboard sections. While this change in chord 
does alter rotor solidity, the reader is reminded that no constraint on rotor solidity has been imposed in the 
current study. However, the 1/4-chord location is invariant to changes in chord resulting from altering blade 
taper. Notice also the upsweep to the trailing edge of the rotor tip at r/R = 1.0. 

Step D shape results are recorded in Figure 6. Again, the modifications are observed over the entire blade 
radius. Chord changes through the taper variable have been allowed without regard to changes in rotor solidity. 



Figure 3. Step A: Composite Weight 2:1, Optimize Thickness, Camber, Collective. Collective perturbation removed and 
vertical scale exaggerated for comparison. 
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Figure 4. Step B: Composite Weight 2:1, Optimize Thickness, Camber, Twist, Collective. Collective perturbation 
removed and vertical scale exaggerated for comparison. 



Figure 5. Step C: Composite Weight 2:1, Optimize Thickness, Camber, Twist, Taper, Collective. Collective perturbation 
removed and vertical scale exaggerated for comparison. 
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r/R=1.00 


Baseline 

Design 


Figure 6. Step D: Composite Weight 1:2, Optimize Thickness, Camber, Twist, Taper, Collective. Collective perturbation 
removed and vertical scale exaggerated for comparison. 


The design history for Steps A-D are represented in Figures 7-10 where the thrust constraint, Ct, was held 
above the baseline value. Here Ct, the objective for each design point FM and 77 , and the composite objective 
function are plotted against design cycle. In each plot, the thrust constraint minimum values are also shown 
for reference. For each design step, both objectives increase rapidly during the first few design cycles. Recall 
that the optimization problem is formulated as a minimization of the composite objective function hence it 
exhibits inverse behavior as it decreases with design cycle. 

VII. Conclusion 

A discrete adjoint-based methodology for performing design optimization of isolated rotor problems that 
appear as steady flows in a noninertial reference frame has been used to drive a multi-point optimization 
of a tilt-rotor configuration. A composite multi-point objective function was constructed using performance 
weighting between the two design points and the optimization was subject to constraints on the thrust at 
each design point. Data was presented for optimization considering a constant thrust and a minimum thrust 
constraint. The optimization successfully increased the hover figure of merit and the airplane-mode propulsive 
efficiency in all cases considered. To better understand the impact of the chosen design variables, a four step 
approach was taken to incrementally increase the geometric degrees of freedom in the design. The results 
illustrated the need to impose a blade loading constraint to limit increases in cruise Ct when designing to 
fixed mission requirements. Nevertheless, the optimization showed further improvement for the TRAM rotor 
previously designed with traditional techniques. 
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Figure 7. Step A: Composite Weight 2:1, Optimize Thickness, Camber, Collective. 
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Figure 8. Step B: Composite Weight 2:1, Optimize Thickness, Camber, Twist, Collective. 
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Step C: Composite Weight 2:1, Optimize Thickness, Camber, Twist, Taper, Collect! 
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Figure 10. Step D: Composite Weight 1:2, Optimize Thickness, Camber, Twist, Taper, Collective. 
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